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Abstract 

The molecular mechanism that maintains the pluripotency of embryonic stem cells (ESCs) is not well understood but may 
be reflected in complex biological networks. However, there have been few studies on the effects of transcriptional and 
post-transcriptional regulation during the development of ESCs from the perspective of computational systems biology. In 
this study, we analyzed the topological properties of the "core" pluripotency transcription factors (TFs) 0CT4, S0X2 and 
NANOG in protein-protein interaction networks (PPINs). Further, we identified synergistic interactions between these TFs 
and microRNAs (miRNAs) in PPINs during ESC development. Results show that there were significant differences in centrality 
characters between TF-targets and non-TF-targets in PPINs. We also found that there was consistent regulation of multiple 
"core" pluripotency TFs. Based on the analysis of shortest path length, we found that the module properties were not only 
within the targets regulated by common or multiple "core" pluripotency TFs but also between the groups of targets 
regulated by different TFs. Finally, we identified synergistic regulation of these TFs and miRNAs. In summary, the synergistic 
effects of "core" pluripotency TFs and miRNAs were analyzed using computational methods in both human and mouse 
PPINs. 
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Introduction 

The capacity to difTerentiate into different cell types, a property 
known as pluripotency, is a defining property of embryonic stem 
cells (ESCs). ESCs are derived from the inner cell mass of the 
mammalian blastocyst [1,2]. Pluripotency may be conferred on 
somatic cells following their fusion with ESCs [3]. During this 
process, the transcription factor (TF) NANOG is specifically 
expressed, and this may facilitate fusion-induced pluripotency [4] . 
Moreover, human and mouse fibroblasts can be reprogrammed 
into ES-like cells which are called induced pluripotent stem cells 
(iPS) by forced expression of other TFs (OCT4, SOX2, Klf4, and 
c-Myc) [5-7] . The quality of iPS is enhanced upon selection of 
cells that express endogenous OCT4 or NANOG [8,9]. Recently, 



Deng et al. reprogramed somatic cells into pluripotent cells using a 
combination of seven smaU-molecule compounds and called them 
CiPS [10]. Epigenetic modifications (DNA methylation, histone 
modification, miRNAs and other methods of epigenetic regula- 
tion) have also been found to play important roles in the 
maintenance of 'stemness' [11-13]. These results indicate that in 
addition to the genetic factors affecting the maintenance of 
pluripotency, complex epigenetic factors are also involved in the 
transformation of ESCs. In order to understand the mechanism by 
which pluripotency is established and maintained in ESCs, further 
effort win be required to research all aspects of the properties of 
molecules and their complex interactions in the biological 
networks which are involved in transcriptional and post-transcrip- 
tional regulation. 
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According to previous studies, a small set of TFs, including 
OCT4, SOX2 and NANOG comprise the "core" pluripotency 
factors in ESCs [14]. OCT4 has long been considered to play 
essential roles in maintaining pluripotency in vivo and in vilro 
[15]. In fact, the concentration of OCT4 is crucial for 
pluripotency: reduced expression evokes trophoectoderm devel- 
opment, whereas enhanced expression leads to primitive endo- 
derm differentiation [16]. As a transcriptional partner of OCT4, 
SOX2 assembles on regulatory elements of target genes together 
with OCT4 to collaborate in transcriptional control, without 
directiy interacting with OCT4 protein [17]. The function of 
NANOG is to promote the self-renewal of ESCs and alleviate the 
requirement for Leukemia Inhibitory Factor (LIE) [18,19]. Among 
OCT4 targets, about half are associated with SOX2. Further- 
more, more than 90% of the target genes shared by OCT4 and 
SOX2 are also associated with NANOG [20]. Based on the above 
results, some researchers have constructed biological networks that 
involve these TFs, and analyzed their properties during the 
development of ESCs [21-24]. In a word, as key factors in the 
maintenance of the pluripotency and self-renewal of ESCs, OCT4, 
SOX2 and NANOG coordinate the regulation of downstream 
genes. 

Besides the traditional genetic impacts on the maintenance of 
ESC pluripotency, epigenetic regulation is also involved in the 
process of ESC development. In particular, more and more studies 
have found that miRNAs play important roles during the 
development of ESCs [25-27]. miRNAs are endogenous single 
strand non-coding RNAs which can inhibit target mRNA 
expression in a post-transcriptional manner [28]. It is character- 
istic of miRNAs that they regulate target genes in a minor manner 
and show temporal and spatial specificity. They may form a 
complex interaction network with other biological molecules in 
vivo. However it is not clear how the target genes of "core" 
pluripotency TFs regulate ESC development synergisticaUy with 
miRNAs. 

In this study, we identified protein-protein interaction networks 
(PPINs) and analyzed the topological properties of the target genes 
of OCT4, SOX2 and NANOG in human and mouse ESCs. 
Further, we explored the effects of miRNAs on the post- 
transcriptional regulation of the target genes of these three "core" 
pluripotency TFs. We found that the centrality of "core" 
pluripotency transcription factor target genes is higher than that 
of randomly selected genes in PPINs. Furthermore, when genes 
are regulated by more "core" pluripotency TFs, they show more 
properties of centrality. The target genes regulated by both 
transcriptional and post-transcriptional methods also have higher 
centrality properties in PPINs. These results indicate that there are 
both the complex interactions between different "core" pluripo- 
tency TFs during ESC development within transcriptional levels 
and the interactions occur across both transcriptional and post- 
transcriptional levels in biological networks. 

Materials and Methods 

Dataset of transcription factor targets 

In order to obtain comprehensive target datasets of "core" 
pluripotency TFs in ESCs, we manually collected related articles 
in PubMed. Finally, 10 articles and their corresponding datasets 
were extracted and used in this study (Table 1). The human 
database contained 3,949 entries, including 623 targets of OCT4, 
1,436 targets of SOX2 and 1,886 targets of NANOG (Figure 1 
left). The mouse database contained 25,222 entries, including 
12,637 targets of OCT4, 5,971 targets of SOX2 and 6,614 targets 



(A) OCT4 SOX2 (B) OCT4 SOX2 




NANOG NANOG 



Figure 1. Venn diagram showing the targets of "core" 
pluripotency TFs. (A) Human targets of TFs. (B) Mouse targets of TFs. 
doi:10.1371/journal.pone.0105180.g001 

of NANOG (Figure 1 right). The detailed targets of these three 
TFs are shown in Table SI and Table S2. 

Datasets of Protein-Protein interactions 

In order to avoid any bias due to the data source, protein- 
protein interaction data were also downloaded from two different 
databases: the Biological General Repository for Interaction 
Datasets (BioGRID) version 3.2.110 (http://theBioGRID.org/) 
and the Human Protein Reference Database (HPRD) (http:// 
www.hprd.org/) [29,30]. We then removed duplicated edges and 
selfloops using Cytoscape and analyzed topological properties 
using the NetworkAnalyzer tools in Cytoscape [31-33]. The 
datasets from BioGRID contained 9,698 nodes with 52,284 edges 
and 4,281 nodes with 7,415 edges (excluding pure high 
throughput experimental data) in human and mouse respectively. 
The HPRD dataset contained 9,453 nodes with 36,867 edges 
(excluding pure high throughput experimental data). 

Dataset of miRNA targets 

Targets of miRNAs were downloaded from two different 
databases, which cover both human and mouse species. The first 
one is the miRecords dataset (http:/ /mirecords. biolead.org/), 
which includes 284 miRNAs, 1,101 targets with 2,087 edges in the 
human and 145 miRNAs, 266 targets with 442 edges in the mouse 
[34]. The second database is TarBase (http://diana.cslab.ece. 
ntua.gr/DianaToolsNew/index.php), which includes 111 miR- 
NAs, 862 targets with 1,093 edges in the human and 44 miRNAs, 
75 targets with 104 edges in the mouse [35]. The targets of 
miRNAs listed in miRecords consist of experimentally verified 
targets and predicted targets which are an integration of predicted 
miRNA targets produced by 1 1 of the following miRNA target 
prediction tools: DIANA-microT, Microlnspector, miRanda, 
MirTarget2, miTarget, NBmiRTar, PicTar, PITA, RNA22, 
RNAhybrid, and TargetScan/TargertScanS. Since TarBase only 
includes experimentally verified miRNA targets, their scale of 
targets is less than those in miRecords. In order to obtain robust 
results, all the datasets used in this work were experimentally 
verified. 

Analysis of protein interaction network topological 
properties 

Analyzing the topological properties of PPINs not only reveals 
the complex molecular interaction pathways, but also provides 
reference points for the detection of important transcriptional 
factors and downstream targets involved in maintaining the 
pluripotency of ESCs. Many topological properties of PPINs were 
used in this analysis. In this study, the Average Shortest Path 
Length (ASPL), Betweenness (BC), Closeness, Clustering Coeffi- 
cient (CC), Degree, Eccentricity, Neighborhood Connectivity 
(NC), Radiality, Stress and Topological Coefficient (TC) were 
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Table 1. Targets of "core" pluripotency TFs in ESCs. 



Species 


Years 


PubMed ID 


OCT4-targets 


SOX2-targets 


NANOG-targets 


human 


2005 


16153702 [20] 


623 


1279 


1687 


human 


2008 


18443585 [46] 


0 


734 


988 


mouse 


2005 


16518401 [47] 


778 


0 


1027 


mouse 


2008 


18555785 [48] 


4369 


2941 


2612 


mouse 


2008 


18358816 [21] 


0 


819 


1284 


mouse 


2008 


18692474 [25] 


4320 


3380 


3114 


mouse 


2008 


18347094 [49] 


2151 


0 


1936 


mouse 


2010 


20362542 [50] 


10 


10 


12 


mouse 


2011 


21477851 [51] 


7518 


0 


0 


mouse 


2013 


23582322 [52] 


4417 


1699 


2492 



doi:l 0.1 371 /journal.pone.01 051 80.t001 



used to analyze the targets of OCT4, SOX2, NANOG and 
miRNAs (TarBases and miRecords) in the PPINs of the BioGRID 
and HPRD databases (Table 2). All the analyses were processed 
with Cytoscape and its NetworkAnalyzer tools. 

All the above topological properties can be used to measure the 
centrality of nodes in biological networks. General speaking, the 
higher the centrality of one node, the more important roles it plays 
in biological networks. For detailed description, we took some 
properties as examples to illustrate their meanings. ASPL is 
defined as the average shortest path length between a node and aU 
the nodes in biological networks. Closeness centrality is defined as 
the reciprocal of the average shortest path length of one node 
which can be used as a measure of how fast information spreads 
from a given node to all other reachable nodes in biological 
networks. In undirected biological networks (such as PPINs), CC 
of a node is defined as the proportion of the observed connections 
between the neighbors of this node against the maximum number 
of possible connections among them. CC is used to indicate the 
close extent of the local neighborhood of one node. Degree is one 



simplest and most used topological index, which is defined as the 
number of nodes directly connected to a given node. TC is a 
relative measurement of the tendency of one node in biological 
networks to have shared interactive partners with other nodes. For 
more in-depth interpretation of these concepts, one can get the 
exact definitions of these topological properties from Table 2. 

Results 

The targets of "core" pluripotency TFs 

A total of 2,512 human pluripotency "core" TF targets were 
identified in the collected articles. Among these targets, 42.9% 
were shared by at least two TFs, corresponding to 1,017 targets. 
The number of OCT4, SOX2 and NANOG targets was 623, 
1,435 and 1,885 respectively, while the proportion shared by other 
TFs was 77.5%, 69.5% and 54.5% respectively (Figure lA). 
Similar results were found in the mouse species. The total number 
of TF targets was 15,714, including 6,393 targets that were shared 
by two or three other TFs, accounting for 40.7% of the total 



Table 2. Definitions of the topological properties. 





Property 


Function 


Description 


Average Shortest Path Length (ASPL) 


ASPL 


The average number of steps for all shortest paths. 


Betweenness Centrality (BC) 


B,= y ^ 


(7k,j denotes shortest paths between node pairs K and j, (Jkij denotes that 
pass through the node 1. 


Closeness Centrality (CC) 


avg{L{n,m)) 


L(n,m) is the length of the shortest path between two nodes n and m. The 
Closeness centrality of each node is a number between 0 and 1. 


Clustering Coefficient (CC) 


^//(d/-l) 


dj is the number of neighbors of i, and e| is the number of connected pairs 
between all neighbors of i. 


Degree 


di 


The number of links to node i. 


Eccentricity 


E 


The maximum node eccentricity (E) can be desicribed as the network 
diameter, that is the largest distance between two nodes. 


Meighborhood Connectivity (NC) 


d 


It is defined as the average connectivity of all neighbors of the node. 


^adiality 


R = D-ASPL+1 


This attribute is a node centrality index computed by the diameter (D) o1 
a node n's the connected component plus 1 and subtracting the average 
shortest path length (ASPL). 


Stress 


s 


s of a node n is the number of shortest paths passing through n. 


Topological Coefficient (TC) 


, _ avg(J{iJ)) 


J(i, j) is the number of neighbors shared between the nodes i and j, plus 
one if there is a direct link between i and j. avg(J(iJ)) is the average value 
of J(i, j). dj is degree of node i. 
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0CT4 S0X2 NANOG 




□ TF-miRXA ■ TF- non-mi RXA □ non-TF-non-mi RNA 



Figure 2. The distribution of targets of "core" pluripotency TFs in liuman BioGRID. (A) miRNA targets obtained from miRecords. (B) 
Targets of miRNAs from TarBase. 
doi:1 0.1 371 /journal.pone.01 051 80.g002 



number. Among the TF targets, 12,636 were OCT4 targets, 5,970 
were SOX2 targets and 6,613 wee NANOG targets. Of these, 
45.0%, 86.8% and 76.2%, respectively, were shared by other TFs 
(Figure IB). As the results show, the target numbers and 
proportions that shared by different TFs were differences between 
human and mouse. The reason may be contributed from the 
research depth on human and mouse. As our manuscript show 
(Table 1), compared with the up to 8 literatures in mouse, there 
were only 2 literatures that contain at least two of the three "core" 
TF targets in human. As one extremely example, there was only 1 
literature that contain OCT4 targets in human. In order to 
overcome the dataset bias and enrich the information in human, 
more works should be done for the genome wide target detection 
of these three "core" TFs in human. 

Mapping targets into PPINs 

We mapped the targets of the three TFs and miRNAs into the 
PPINs and obtained sub-networks consisting only of "core" 
pluripotency TF targets and their direct neighbors (Figure SI to 
S9). From both BioGRID and HPRD results, it was evident that 
the proportions of TF-miRNA targets (targets which are regulated 
by both these three "core" TFs and miRNAs) and TF-non- 
miRNA targets (targets which are regulated only by these three 
"core" TFs but not regulated by miRNAs) were smaller compared 
with the proportion of non-TF-non-miRNA genes (genes which 
are neither regulated by these three "core" TFs nor regulated by 
miRNAs) in both human and mouse (Figure 2, 3 and 4). 



Analysis of topological properties 

The topological properties were analyzed with Cytoscape 
version 3.0.2, especially the NetworkAnalyzer tools. To begin 
with, the most connected components of protein-protein interac- 
tion networks were extracted from BioGRID, including 9,552 
nodes with 52,202 edges in the human and 3,831 nodes with 7,123 
edges in the mouse. The same analysis identified 9,205 nodes and 
36,748 edges in the most connected component of the HPRD 
database. Finally, topological properties were analyzed with 
NetworkAnalyzer and filtered using strict statistical parameters 
[P values<0.05 with t test, processed with R version 3.0.2). 

Comparison of centrality properties between TF-targets 
and non-TF-targets 

Following, we compared the topological properties between TF- 
targets (genes that are regulated by the "core" TFs, including 
OCT4, SOX2 and NANOG) and non-TF-targets (genes that are 
not regulated by any of the "core" TFs, including OCT4, SOX2 
and NANOG). First we analyzed human PPINs. The results 
showed that the ASPL of SOX2 -targets and NANOG-targets was 
shorter compared with non-targets, in both the BioGRID and 
HPRD datasets. We also found that radiahty of SOX2 -targets and 
NANOG-targets was greater than non-targets in both BioGRID 
and HPRD resources. Furthermore, the degree of the NANOG- 
targets was also significantly different compared with non- 
NANOG-targets in BioGRID and HPRD, indicating that many 
proteins are connected with NANOG-targets. Without consistent 
significant results, the SOX2-targets were only found to have 
higher degree values compared with non-SOX2-targets in HPRD, 
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0CT4 S0X2 NANOG 




□ TF-miRXA ■ TF- non-mi R\A □ non-TF-non-mi R\A 



Figure 3. The distribution of targets of "core" pluripotency TFs in HPRD. (A) miRNA targets obtained from miRecords. (B) Targets of miRNAs 
from TarBase. 

doi:10.1371/journal.pone.0105180.g003 




□ TF-miR\A ■ TF- non-mi RXA □ non-TF-non-mi RXA 



Figure 4. Tlie distribution of targets of "core" pluripotency TFs in mouse BioGRID. (A) miRNA targets obtained from miRecords. (B) 
Targets of miRNAs from TarBase. 
doi:1 0.1 371 /journal.pone.01 05 1 80.g004 
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while similar results were not found in BioGRID. With regard to 
OCT4-targets, BC difTered significantly between OCT4-targets 
and non-OCT4-targets in both BioGRID and HPRD databases, 
indicating that the shortest paths going through OCT4 targets 
were more than a random choice. This indicates that OCT4 
targets may be internal module proteins and are more likely to 
locate in the hub position in networks. In summary, a certain 
degree of higher centrality in PPINs was found in human "core" 
pluripotency targets compared with non-TF-targets (Table 3). 

Similar results were also obtained in the mouse. For NANOG- 
targets, 6 measurements were found to differ significantly from 
those in non-NANOG-targets, including ASPL, Closeness, 
Degree, NC, Radiality and Stress. For SOX2-targets, 5 measure- 
ments in total were significantly different compared with non- 
SOX2-targets: ASPL, Closeness, Degree, NC and Radiality. For 
OCT4-targets, 5 measurements were found to differ from noii- 
OCT4-targets, including BC, Degree, NC, Stress and TC. Taking 
these measurement results together, the target genes of "core" 
pluripotency transcription factors show higher centrality properties 
in mouse PPINs (Table 4). 

Consistency analysis of multiple "core" pluripotency IF 
regulations 

Through the analysis of the distributions of "core" pluripotency 
TF targets, we identified many genes regulated by at least two TFs 
(Figure 1). This result indicates that these TFs may be involved in 
complex interactions and execute similar functions synergistically 
as cells progress along the pathway of ESC development. To 
investigate this further, we continued to explore cooperation 
between the TF regulators through their topological properties. As 
expected, no difference in the centrality properties was found 
between ITF-targtes, 2TF-targets and 3TF-targets in human 
PPINs, whichever protein-protein interaction datasets were used 
(Table 5). 

Similar results were also found in the mouse. Among the ten 
centrality properties, none of them was found to differ between 
ITF-targets and 2TF-targets. When compared with 3TF-targets, 
only ASPL, Closeness and Radiality were found to be different 
from those of 2TF-targets. The greatest diversity was found 
between the groups of the 1 TF-targets and 3TF-targets, where we 
found differences in five measurements: ASPL, Closeness, 
Eccentricity, Radiality and Stress (Table 6). 

Modularity within the inner "core" pluripotency TFs with 
shortest path length analysis 

When we investigated the regulation of multiply "core" 
pluripotency TFs during the development of ESCs in PPINs, 
one question was triggered about whether there are closer 
relationships between targets of these TFs. One hypothesis 
suggests that the connections within and between TF targets 
should be closer than between other genes in PPINs. In other 
words, module properties are expected within and between the 
"core" pluripotency TFs. In order to verify this hypothesis, we 
performed shortest path length (SPL) analysis across human and 
mouse species. We found that the background averages of PPIN 
SPLs were 4.227, 4.201 and 5.125 in HPRD, human BioGRID 
and mouse BioGRID respectively. First, the smaller SPLs were 
detected within TF targets compared with the background SPLs. 
As shown in Table 7, 8 and 9, the SPLs of TF targets, including 
OCT4, SOX2 and NANOG, were aU smaUer than those of other 
proteins in the PPINs, no matter which source of PPIN data was 
used. Second, the common targets of at least two TFs showed 
smaUer SPLs compared with other proteins in the PPINs (Table 7, 
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Figure 5. Analysis of topological properties of OCT4-miRNA-targets in human PPINs. (A) Comparison of 0CT4-miRNA and OCT4-non- 
miRNA targets. (B) Comparison of 0CT4-miRNA and non-0CT4-non-miRNA targets. 
doi:1 0.1 371 /journal.pone.01 051 80.g005 
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Figure 6. Analysis of topological properties of SOX2-miRNA-targets in human PPINs. (A) Comparison of S0X2-miRNA and S0X2-non- 
miRNA targets. (B) Comparison of S0X2-miRNA and non-SOX2-non-miRNA targets. 
doi:1 0.1 371 /journal.pone.01 051 80.g006 
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Figure 7. Analysis of topological properties of NANOG-miRNA-targets in human PPINs. (A) Comparison of NANOG-miRNA and NANOG- 
non-miRNA targets. (B) Comparison of NANOG-miRNA and non-NANOG-non-miRNA targets. 
doi:1 0.1 371 /journal.pone.01 051 80.g007 



8 and 9). When the number of TFs in combination was 2, the 
corresponding P values of the t tests were 4.99E-5, 0.001 and 
3.90E-286; while when the number was 3, the P value was 7.98E- 
258 for values from the HPRD database (Table 7). Similar results 
were found for other databases (Table 8 and 9). Finally, we 
compared the distances between groups of different TF targets 
with SPLs and found that they were significantly different, 
especially when the group comprised the common targets of three 
TFs (Table 7, 8 and 9). From the HPRD and mouse BioGRID 
databases we found that the SPLs of three TFs were smaller than 
those of most other groups, but similar results were not observed in 
the PPIN of the human BioGRID. Summarizing the above results, 
we found module properties not only within the targets regulated 
by common or multiple "core" pluripotency TFs but also between 
the groups of targets regulated by different TFs. 

Post-transcriptional regulation effects on "core" 
pluripotency TF targets 

Based on the evidence that epigenetic regulation may play 
important roles in ESC development, we attempted to analyze the 
effects of post-transcriptional regulation on the targets of "core" 
pluripotency TFs in PPINs. As typical forms of post-transcriptional 
regulation, miRNA targets were obtained for further analysis from 
two different sources, including miRecords and TarBase. Results 
show that there were some significant different characteristics 
between TF-miRNA targets and TF-non-miRNA targets in the 
both two PPINs from human BioGRID and HPRD databases. 
The differences were found in both miRecords and TarBase. First, 
we found that OCT4-miRNA targets were different from OCT4- 
non-miRNA targets in several ceiitrality properties, including 



Degree, Stress and TC (Figure 5). Second, observation of both 
SOX2-miRNA targets and NANOG-miRNA targets revealed that 
they differed from TF-non-miRNA targets in several measure- 
ments. The measurements that differed between SOX2-miRNA 
and SOX2-non-miRNA targets were ASPL, Closeness, Degree, 
Eccentricity, Radiality, Stress and TC (Figure 6). And the 
measurements, including ASPL, Clossiiess, Degree, Eccentricity, 
Radiality, Stress and TC, were different between NANOG- 
miRNA targets and NANOG-non-miRNA targets (Figure 7). 
Furthermore, these results were observed across both the 
BioGRID and HPRD protein-protein interaction databases and 
the miRNA target databases miRecords, TarBase. 

Further, in order to identify the more obvious effects of miRNA 
regulation on the development of ESCs, we compared TF-miRNA 
target genes to non-TF-non-miRNA target genes in the human. 
Interestingly, as expected we found a significant difference 
between TF-miRNA and non-TF-non-miRNA targets on the 
topological properties of components of PPINs. First, OCT4- 
miRNA targets were found to have higher centrality properties 
compared with non-OCT4-non-miRNA genes in PPINs, with 
higher values of ASPL, BC, Degree, Stress, Radiality and TC 
(Figure 5). Second, the measurements which were higher in 
SOX2-miRNA targets compared with non-SOX2-non-miRNA 
genes were ASPL, BC, Closeness, Degree, Eccentricity, Radiality, 
Stress and TC (Figure 6). Third, up to 9 measurements differed 
between NANOG-miRNA targets and non-NANOG-non- 
miRNA genes in PPINs, excluding CC (Figure 7). These results 
were supported by all the different data sources, including 
BioGRID, HPRD, miRecords and TarBase. 
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Figure 8. Analysis of topological properties of OCT4-miRNA-targets in mouse PPIN. (A) Comparison of 0CT4-miRNA and 0CT4-non-miRNA 
targets. (B) Comparison of 0CT4-miRNA and non-OCT4-non-miRNA targets. 
doi:1 0.1 371/journal.pone.01 051 80.g008 



In order to overcome species bias, we conducted the same 
experiments in mouse ESCs. Interestingly, the results showed the 
same tendency as above in data from all the dataset sources. 
Firstly, the TF-miRNA targets also showed higher centrality 
properties compared with TF-non-miRNA targets (Figure 8, 9 
and 10). The numbers of measurements that were different 
between TF-miRNA targets and TF-non-miRNA targets were up 
to 7, 9 and 5 in OCT4, SOX2 and NANOG respectively. The 
centrality properties that were not found to differ between OCT4- 
miRNA and OCT4-non-miRNA targets were CC, NC and TC. 
For SOX2, only NC showed no difference. The higher 
measurements in NANOG-miRNA targets included ASPL, 
Closeness, Degree, Eccentricity and Radiality. Second, the TF- 
miRNA targets differed from the non-TF-non-miRNA target 
genes in PPINs (Figure 8, 9 and 10). The numbers of measure- 
ments that were different between TF-miRNA targets and non- 
TF-non-miRNA targets were up to 9, 10 and 7 in OCT4, SOX2 
and NANOG respectively. The measurements that were not found 
to differ were CC and BC, CC and TC in OCT4 and NANOG 
respectively. In summary, miRNAs play important roles during 
the development of ESCs and participate in complex interactions 
with "core" pluripotency TFs from a biological system perspective. 

Discussion 

In this study, we globally analyzed the topological properties of 
targets of the "core" pluripotency TFs in PPINs, including OCT4, 
SOX2 and NANOG. In addition, the post-transciptional effects of 
miRNAs on these TFs were also analyzed in both human and 
mouse. Up to ten topological properties were included in the 



analysis process, including Shortest Path Length, Betweeness, 
Closeness, Cluster Coefficient, Degree, Eccentricity, Neighbor- 
hood Connectivity, Radiality, Stress and Topological Coefficient. 
AH the above analyses were processed in three protein-protein 
interaction datasets (HPRD, human BioGRID and mouse 
BioGRID), two miRNA target datasets (miRecords and TarBase) 
and two species (human and mouse). Though there were different 
dataset scales of the three "core" TF targets between different 
databases and species, the common characteristic of these three 
"core" TFs in biological networks were still observed in our 
research. The use of several data sources and measurements in this 
study ensures the robust nature of the results obtained. Besides, all 
the above analysis was based on ESCs datasets. Consider the 
similar property of "stemness" with ESCs in many other types of 
stem cells, especially induced pluripotent stem cells, we infer that 
these three "core" TFs will have similar roles and characteristics in 
biological networks. With the increasing of high throughput 
datasets about targets of the "core" pluripotency TFs in other 
types of stem cells, similar research should be processed and 
compared with current results. 

We found significant differences in centrality properties between 
"core" pluripotency TF-targets and non-TF-targets in PPINs. 
These results were widespread in HPRD, human BioGRID and 
mouse BioGRID datasets. The numbers of centrality properties 
were 6 and 8 in human and mouse respectively. The former 
contained ASPL, BC, Degree, Eccentricity, Radiality and Stress, 
while the latter comprised ASPL, BC, Closeness, Degree, NC, 
Radiality, Stress and TC. Comparing the two results, we found 
that ASPL, BC, Degree, Radiality and Stress were robust and only 
CC among the 10 measurements could not be detected, which is 
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Figure 9. Analysis of topological properties of SOX2-miRNA-targets in mouse PPIN. (Aj Comparison of S0X2-miRNA and S0X2-non-miRNA 
targets. (B) Comparison of S0X2-miRNA and non-S0X2-non-miRNA targets. 
doi:1 0.1 371/journal.pone.01 051 80.g009 



used to judge the close link of node neighborhoods in biological 
networks. The reason why CC did not appear in the analysis was 
not clear. It may be because the targets of "core" pluripotency TFs 
perform their functions in relative isolation which may help to 
avoid harm to the complex environment in vivo during the 
development of ESCs [36,37]. With the higher central properties, 
these results indicate that targets of these three "core" TFs play 
more important roles than random genes during the development 
of ESCs. 

We found synergistic regulation of multiple "core" pluripotency 
TFs during the development of ESCs. As we hypothesized, no 
difference in topological properties was found between ITF- 
targets, 2TF-targets and 3TF-targets in human PPINs, including 
HPRD and human BioGRID. The same tendency was also found 
in mouse PPINs from BioGRID. Based on these results, it can be 
seen that although the number of TFs regulating common target 
genes increases from 1 to 3, their centrality properties do not 
increase accordingly. This indicates that pluripotency related 
genes may be regulated by 1 or 2 or even 3 TFs, but the genes are 
no different from the biological system viewpoint. In other words, 
the number of TFs that regulate the common pluripotency genes is 
not important. This means that the seemingly unnecessary TFs 
may provide compensatory regulation in the molecular pathways 
of ESC development [38]. Such compensatory regulation wiU 
conversely enhance the status of the common targets in a 
biological process. Through the consistent regulation of these 
TFs, the maintenance of the pluripotency of ESCs is rendered 
more reliable. As an example, the synergistic regulation of histone 
deacetylase I (HDACl) by OCT4, SOX2 and NANOG plays 



important roles not only in the development of ESCs but also in 
the growth of tumor cells [39-41]. 

The module properties within the inner "core" pluripotency 
TFs were detected through the analysis of shortest path length in 
PPINs. Through these experiments, we found that the shortest 
path length of targets regulated by common TFs were smaller than 
those of randomly selected background values. This indicates that 
the genes regulated by common "core" pluripotency TFs are in 
close communication in PPINs and that this will be helpful in 
synergisticaUy and quickly maintaining the pluripotency of ESCs 
[42] . Next, the shortest path lengths of target genes regulated by 
multiple TFs were further analyzed. Results show that these genes 
also have smaller shortest path lengths between each other. Third, 
we took the analysis of the extent of contact between groups of 
different TF targets and found that they were different, especially 
when the group comprised the common targets regulated by three 
TFs. In summary, we identified module properties not only within 
the targets regulated by common or multiple "core" pluripotency 
TFs but also between the groups of targets regulated by different 
TFs. The genes within module may have close and quickly 
information flowing which will help them to implement the 
common function of maintaining pluripotency during the devel- 
opment of ESCs. 

We also found that miRNAs play important roles in the 
regulation of "core" pluripotency TFs during the development of 
ESCs as a way of post-transcriptional regulation in PPINs. The 
diflference of centrality properties observed between TF-miRNA 
targets and TF-non-miRNA targets was found in both human and 
mouse PPINs. Further, the differences in topological properties 
between TF-miRNA targets and non-TF-non-miRNA target 
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Figure 10. Analysis of topological properties of NANOG-miRNA-targets in mouse PPIN. (Aj Comparison of NANOG-miRNA and NANOG- 
non-miRNA targets. (B) Comparison of NANOG-miRNA and non-NANOG-non-miRNA targets. 
doi:10.1371/journal.pone.0105180.g010 



genes were more obvious in PPINs. These different centrality 
properties show difTerent close extent and different importance in 
computational biological network view. Considering the effects 
that miRNAs impose to the "core" TFs on these properties in 
PPINs, it is observed that there is synergistic regulation between 
themselves. The synergistic regulation of "core" pluripotency TFs 
and miRNAs will enhance the function of targets. In consideration 
of the drastic effects of TFs and the minor regulation of miRNAs, 
the synergism of these three "core" TFs and miRNAs may help to 
achieve their aims of regulation about pluripotency, like one 
machine which has many gear wheel of different size. As an 
example, myocyte enhancer factor2 (MEF2) which is a target of 
both NANOG and miRNA, is an important transcription factor 
regulating the survival and development of many types of cells 



[43,44]. One of its prominent functions is the control of gene 
transcription in cell differentiation. AH of the other genes that 
regulated by NANOG and miRNAs (including Mxil, AridSb, Kit, 
Hoxall, Hoxa7, Mef2c, Gjal, Myc, Hdac4 and Tmsb4x) are 
known to be related with stem cells. To further verify our results, 
we performed analysis of function categories with these ten targets 
in PIR and found that they are development proteins and related 
with transcription regulation (Table 10). Our results clearly reveal 
the effects of epigenetic regulations on the development of ESCs in 
biological networks, findings which are consistent with previous 
studies performed using molecular and cell technology [25,45]. 
This finding may provide candidate pathway for deep detection 
about ESCs molecular mechanism from post transcriptional level. 



Table 10. Function categories analysis about targets of NANOG and miRNAs in PIR. 




Term 


P Value 


Transcription regulation 


8.59E-05 


Transcription 


9.74E-05 


Dna-binding 


7.75E-04 


Proto-oncogene 


0.004849 


Nucleus 


0.005432 


DNA binding 


0.010331 


Developmental protein 


0.048812 


doi:10.1371/journal.pone.0105180.t010 
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